Astronomy & Astrophysics manuscript no. h4435 


February 2, 2008 


(DOI: will be inserted by hand later) 





Pulse-to-pulse intensity modulation and drifting subpulses in 

recycled pulsars 

R. T. Edwards 1 and B. W. S tappers 2 ' 1 



m 
O 
O 



in 



> 

(N 

in 
o 
m 
o 

Or 
i 1 

o ■ 
a: 

C/3 . 



13 



1 Astronomical Institute "Anton Pannekoek", University of Amsterdam, Kruislaan 403, 1098 SJ Amsterdam, The Netherlands 

2 Stichting ASTRON, Postbus 2, 7990 AA Dwingeloo, The Netherlands 

the date of receipt and acceptance should be inserted later 

Abstract. We report the detection of pulse-to-pulse periodic intensity modulations, in observations of recycled pulsars. Even 
though the detection of individual pulses was generally not possible due to their low flux density and short duration, through 
the accumulation of statistics over sequences of 10 s — 10 6 pulses we were able to determine the presence and properties of 
the pulse-to-pulse intensity variations of six pulsars. In most cases we found that the modulation included a weak, broadly 
quasi-periodic component. For two pulsars the sensitivity was high enough to ascertain that the modulation phase apparently 
varies systematically across the profile, indicating that the modulation appears as drifting subpulses. We detected brighter than 
average individual pulses in several pulsars, with energies up to 2-7 times higher than the mean, similar to results from normal 
pulsars. We were sensitive to giant pulses of a rate of occurrence equal to (and in many instances much lower than) that of PSR 
B1937+21 at 1400 MHz (~ 30 times lower than at 430 MHz), but none were detected, indicating that the phenomenon is rare 
in recycled pulsars. 
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1. Introduction 

Pulsar science is perhaps unusual in astronomy, in that progress 
in the theoretical understanding of the emission mechanism is 
not limited by a lack of data, but rather by an overabundance of 
complex observable pheonomenology. Pulsar signals are seen 
to vary on every ob servable timesc ale from nanoseconds (e.g . 
lHankins et alJ2003l) to decades (e.g. lWeisberg & Tavlorl2 002). 
with variations associated with such phenomena as microstruc- 
ture, subpulses, periodic/drifting subpulses, giant pulses, pulse 
nulling, mode changing, interstellar scintillation and geodetic 
spin precession. All of these phenomena have been important 
in shaping ideas concerning the radio pulsar emission process, 
however to date a complete explanation of radio pulsar emis- 
sion remains elusive. 

This discovery of PSR B1937+21 (Bac ker et all 1 1982I) 
posed further challenges due to its very short spin period of 
1.56 ms. As more "millisecond pulsars" (MSPs) were discov- 
ered, the rapid spin and low inferred surface magnetic field 
strength grew to be understood as the end product of accre- 
tion from a binary companion, "rec ycling" pulsars that had ear - 
lier evolved into radio silence (e.g. Phinne v~& Kulkarnill 994). 
However, the question of how the radio emission mechanism 
is related to that of ordinary pulsars remained open. With a 
difference of 2-3 orders of magnitude in rotation period and 
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3-4 orders of magnitude in magnetic field strength, a common 
emission mechanism might not be expected, yet a good argu- 
ment can be made for this through the comparison of (polari- 
metric) average pulse profile morphologies and their frequency 
dependence ( Kr amer & X ilouris 2000 and references therein). 

As noted above, the true pulsar signal varies on all 
timescales and much has been learned about "ordinary" 
pulsars from studying aspects of the signal other than 
the basic average profile. This has not been the case 
with MSPs because their typically lower flux densities 
and shorter pulses make the detection of individual pulses 
difficult. An exception has been t he study of so -called 
"giant" pulse s in PSR B 1937+2 1 JWolszczan et all 1 1984 
Bac keil Il995l ICognard et all Il996l iKinkhabwala & Thorsettl 
20001 IVivekanandll2002t IPopov & Steppers! 120031) . and PSR 
B1821-24 jRomani& Johnston! l200ll) . Studies of the prop- 
erties of pulses of normal intensity have thu s far been lim- 
ited to three pulsars. Sallmen & Backer ( 1995) used the 305-m 
Arecibo dish to obtain sensitivity to individual pulses of PSR 
B 1534+ 12 of average and above-average energy and found that 
their distribution in energy was similar to that of some ordi- 
nary pulsars. On the ba sis of the d etection o f pulses of above- 
avera ge intensity, both ljenet et alJ Jl998l) and lVivekanand et alJ 
( 1998) found that the properties of the brightest known MSP, 
PSR J0437-4715, were similar to those of ordinary pulsars, al- 
though to good significance neither nulling nor the presence of 
preferred timescales within individual pulses (e.g. microstruc- 
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ture, subpulses) was detected. Both studies report the detection 
of quasi-periodic pulse-t o-pulse intensity modulat ions with a 
period of ~ 4 pulses, and lVivekanand et al.l (fl998) report (we 
believe potentially erroneously; see Sect. 0} that the modula- 
tion is not associated with drifting subpulses, as it often is in 
ordinary pulsars. Finally. Ijenet et al ] d200 111 used data with an 
average single-pulse signal-to-noise (S/N) ratio less than one 
to show that on a statistical basis and with the exception of gi- 
ant pulses, the emission of PSR B 1937 +21 is extremely stable, 
showing to high significance a complete lack of pulse-to-pulse 
variations and no difference between the shape of individual 
pulses and that of th e average p rofile. 

As shown by Ijenet et all (bOOll) . poor single-pulse S/N 
does not eliminate the possibility of obtaining useful informa- 
tion regarding single pulses. The average profile can of course 
be measured in such cases, but there is no reason why col- 
lection of statistics cannot be extended beyond the first mo- 
ment to include higher-order moments and correlations. Of pri- 
mary interest are the second moment 1 (giving access to the 
modulation index), and the second order correlatio ns. A vari- 
ety of useful correlation statistics can be derived; Ijenet et all 
integrated the single-pulse autocorrelation function 
(ACF) to show that every pulse appeared to be identical to 
the average profile. We extend this approach here to include 
the L ongitude-Resolved Fluctuation Spectrum (LRFS; B acker! 
1970l). the Two-Dimensio nal Fluctuation Spectrum (2DFS; 
Edwards & Stappe rsl2002l) and the Longitud e-Resolved Cross- 
Correlation Function (LRCCF; Popov] Il986l) . Applying these 
techniques to archival observations made at the Westerbork 
Synthesis Radio Telescope (WSRT), we have detected and 
characterised pulse modulation behaviour in several recycled 
pulsars. 

2. Observations and Methods of Analysis 

2.1. Observations 

The data used in this project were selected from archival data 
taken at WSRT. For many pulsars numerous observations were 
available, allowing us to take advantage of fortuitous scintilla- 
tion conditions for substantial enhancements in sensitivity. For 
all observations, dual linear polarization signals from fourteen 
25-m dishes were added using previously determined phase 
and gain factors and the resultant signal was digitally processed 
by the PuMa pulsar backend to form a two-dimensional array 
of total power samples , as a funct i on of time and radio fre- 
quency; for details see Vout e et alJ d2002). In offline process- 
ing, frequency channels containing periodic interference were 
flagged and a de-dispersed time series was produced by com- 
bining remaining channels, aligned using previously published 
dispersion measures. In Table^we list parameters of the obser- 
vations used, including the date of observation as a Modified 
Julian Day (MJD), the observation duration (7\,b s ) in seconds, 
the centre frequency (v) and total bandwidth (Av) in megahertz, 
the number of frequency channels used, the resultant dispersion 

1 Since they worked from voltage samples instead of intensity sam- 
ples, |]enet_et al. (2001) referred to this as the fourth moment. 



smearing (T smea r, near the center frequency) and output sam- 
ple interval (r samp ) in microseconds, and the number of pulses 
recorded . The rightmost four columns refer to the sensitivity, 
see Sect. fJH] 

2.2. Production of a Longitude/Time-dependent Array 

For many techniques of single pulse analysis including those 
used in this work, the data need to be treated as a sequence 
of sampled pulses, with a consistent sampling lattice within 
each pulse. The simplest approach to producing a such two- 
dimensional (i.e. longitude- and time-resolved) representation 
of the pulse sequence is to bin each sample according to accord- 
ing its pulse phase, as calculted using an ephemeris (for which 
purpose we used the TEMPO software package 2 ). When the 
time resolution is close to the interval corresponding to the de- 
sired longitude resolution, a problem arises with this method. 
In the simplest case, where the longitude bin width is chosen 
to correspond to one sample interval (at some point during the 
observation), each pulse as it appears in the binned sequence 
is effectively shifted by an amount corresponding to the offset 
between each sample and the bin centers. This offset advances 
(modulo one bin) by a constant amount each pulse due to the 
fact that the apparent pulse period does not (in general) equal 
an integer number of sample intervals. The pulse shape distor- 
tion (longitude shift) is thus periodic from pulse to pulse, and 
gives rise to a sequence of harmonics in the LRF S. This prob- 
lem was encountered bv lVivekanand et all Jl998l) in their anal- 
ysis of data from PSR J0437-1715, but no attempt was made 
to remove it. 

In this work we avoided the effect by compensating for 
the shift in each pulse using a frequency-domain time shift. 
Formally, we may view the time series as a sampled signal 
that represents the source intensity signal, convolved with a 
function representing the applied integration. By using a se- 
quence of samples that is offset from the required set of bin 
centers (longitude samples) by some fraction of a sample e e 
[-0.5,0.5], and assigning those samples to the nearest bin, 
we arrive at a version of the pulse signal that is apparently 
longitude-shifted by e samples. This can be simply corrected 
by convolution with an offset sine function, which by virtue of 
Fast Fourier Transform algorithms is efficiently performed by 
multiplying the Discrete Fourier Transform (DFT) of the pulse 
by A(vi) = exp(27n'evi) and taking the inverse DFT 3 . If this 
correction is made to each pulse, the artifact is removed com- 
pletely from the signal and all fluctuation statistics can gener- 
ally be interpreted correctly. 

In some circumstances, extension of the approach to ef- 
fect a change in sampling interval (by dropping coefficients 



http://pulsar.princeton.edu/tempo/ 
3 Since this performs a cyclical convolution, it is necessary to dis- 
card some samples from the ends of the sequence. This is most not- 
icable in the presence of strong interference, when the apparent step 
between the last sample and the first can cause ringing at the Nyquist 
freqency when the sampling lattice is offset by this process. We found 
that discarding 32 samples from each end of the result was sufficient 
to avoid this effect. 
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Table 1. Observation parameters and sensitivity 
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Date 
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Av 
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(PSR) 


(MJD) 


(s) 


(MHz) 


(MHz) 




(yUS) 


QiS) 


(10 3 ) 








(< E >) 


J0034-0534 


52216.1 


1800 


328 


10 


512 


63.2 


102.4 


775 


2.0 


1.7 


0.11 


20 


J02 18+4232 


52462.4 


1800 


328 


10 


1024 


140.7 


102.4 


959 


3.6 


3.0 


0.19 


26 


J0613-0200 


52190.1 


9000 


1380 


80 


512 


19.1 


51.2 


2914 


3.1 


2.6 


0.16 


19 


J1012+5307 


52145.7 


1800 


840 


80 


512 


19.7 


102.4 


342 


0.37 


0.31 


0.019 


3.3 




52392.8 


1800 


1380 


80 


512 


4.5 


51.2 


342 


0.2 


0.17 


0.011 


0.78 


J1022+1001 


52004.8 


7200 


328 


10 


512 


47.0 


102.4 


433 


0.29 


0.24 


0.015 


1.5 




52004.6 


7200 


840 


80 


64 


22.4 


51.2 


433 


0.55 


0.47 


0.029 


2.2 




52086.6 


3600 


1380 


80 


512 


5.1 


102.4 


433 


0.21 


0.17 


0.011 


0.62 


J1518+4904 


52101.7 


3600 


1390 


60 


384 


5.6 


204.8 


85 


0.096 


0.081 


0.0057 2 


0.65 


J 1643- 1224 


52242.5 


1800 


1380 


80 


512 


30.8 


102.4 


368 


1.0 


0.85 


0.053 


12 


J17 13+0747 


52544.7 


7200 


840 


80 


512 


35.0 


51.2 


1558 


0.54 


0.46 


0.028 


8.7 




52515.7 


1800 


1190 


80 


512 


12.3 


51.2 


376 


0.12 


0.10 


0.0064 


1.8 




52537.7 


3600 


1700 


40 


256 


4.2 


51.2 


770 


0.11 


0.092 


0.0058 


2.1 




52572.5 


3600 


2240 


80 


512 


1.8 


51.2 


770 


0.27 


0.23 


0.014 


5.2 


J 19 18-0642 


52560.8 


1800 


1380 


80 


512 


13.1 


102.4 


235 


0.58 


0.49 


0.030 


5.8 


B1937+21 


52403.3 


900 


1380 


80 


512 


35.0 


25.6 


576 


0.31 


0.26 


0.016 


5.3 


J2145-0750 


52223.7 


1500 


860 


40 


256 


19.7 


102.4 


91 


0.38 


0.32 


0.020 


1.2 




52232.9 


3600 


1380 


80 


512 


4.4 


102.4 


222 


0.21 


0.18 


0.011 


0.92 



1 :M max = 24576 2 : M lmx = 81920 



in the frequency domain) might be indicated. This is often 
necessary in the analysis of simultaneous multi-observatory 
data sets due to the employment of different sampling inter- 
vals, and may also become important in future high-resolution 
studies of close-orbit binary pulsars where the apparent period 
can change significantly over the course of an observation. We 
suggest the technique be investigated as a means for avoiding 
the deleterious effects o f time-domain re-sampling noticed by 
iKarasterdou et all J20oiR 

In theory, the above procedure allows absolute alignment 
in longitude of observations at different epochs and frequen- 
cies. In practice, the published dispersion measure may derive 
from arbitrary profile alignment procedures, and previously 
published timing ephemerides may not extrapolate beyond the 
date range from which they derive accurately enough for phase 
alignment between different epochs. For these reasons, we did 
not attempt absolute alignment, and our longitude axes include 
an arbitrary offset. 

2.3. Accumulation of Statistics 

The first step of processing after forming the longitude/time ar- 
ray was to correct for an absolute offset due to the system noise 
power. To reduce the effect of slow variations in the system 
temperature we subtracted a running baseline, computed as the 
mean over the surrounding ± 1000 pulses of all samples in a 
defined "off-pulse" longitude interval. Denoting the result as 
S ij where i and j are indices in pulse longitude and pulse num- 



4 Although we also note that for some purposes, including appli- 
cation of the LRCCF between observatories as in K arastergiou et alj 
1200 ll) . the synthesis of synchronous samples is unneccessary and 
should be avoided. 



ber, we then computed the (normalized) average pulse profile 

as 

1 N ~ l 

where N is the number of pulses analyzed. For convenience, at 
this point we re-normalized S u so that the peak bin of /i, had a 
value of 1 . We then computed a version of the data from which 
the average profile was subtracted: 

•V; ; -.V ..;-/';• (2) 

The remaining analysis involved the integration of three 
kinds of correlation statistics that differ in their treatment of the 
longitude dimension. These are the full two-dimensional auto- 
correlation function (2DACF), the longitude-resolved autocor- 
relation function (LRACF) and the longitude-resolved cross- 
correlation function (LRCCF). 

To examine the autocorrelation statistics of the signal while 
keeping computational time and space requirements manage- 
able, we divided the input pulse series into N/M M-pulse 
blocks of data and obtained estimates of the fluctuation power 
spectra (LRFS and 2DFS) by averaging the results of N/M in- 
dividual spectra computed from the squared modulus of the ap- 
propriate DFTs. By application of the Wiener-Khinchin the- 
orem, the corresponding autocorrelation functions were ob- 
tained via the DFTs of the power spectra. The choice of M 
affected our sensitivity to nearly coherent modulation: since be- 
yond this length we summed power spectra, signals with decor- 
relation time-scales longer than M pulses would have been 
detected with reduced significance. Signals with decorrelation 
timescales shorter than M could also have suffered since in that 
case the response would be spread over several coefficients, 
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however here the sensitivity can be recovered by smoothing 
the power spectrum. Our approach was therefore to choose the 
largest M that was practical given the computational and mem- 
ory requirements. 

In computation of the LRFS, each one-dimensional spec- 
trum was obtained using the squared modulus of the coeffi- 
cients of DFTs performed along a constant-longitude column 
of the longitude/time array: 



N/M-l 



LRFS™ 



(i+l)Af-l 



jf Z Z s « e 

k=o i=m 



■2njlyp\/M 



(3) 



To remove the effect of the variance of the system noise, includ- 
ing non-white components due to interference and system tem- 
perature variations not removed in the pre-processing stages, 
we subtracted from each element of the spectrum the corre- 
sponding value from the mean fluctuation spectrum of all des- 
ignated "off-pulse" longitude bins, dropping the superscript of 
Eq.[3]to denote the noise-corrected LRFS. 

We then performed DFTs along constant-longitude 
columns of the LRFS to obtain the LRACF. To avoid perform- 
ing a cyclic convolution, the individual LRFS were actually 
computed over 2M points, of which the elements of one con- 
tiguous half were set to zero. The resulting LRFS elements are 
not independent, so to facilitate proper evaluation of the signifi- 
cance of spectral features, after using it to form the LRACF we 
discarded the odd-numbered coefficients of the LRFS (which 
has the same result as an M-point transform of the non-padded 
data). To correct the LRACF result for edge normalization ef- 
fects we divided the result by 1 - j/M, where is j the lag num- 
ber. The result is: 



LRACF,, = 



1 



N(l - j/M) 



N/M-l (Jc+V)M-l-j 
Z Z S 'i' S 'i(i+j)- 



(4) 



fe=0 



I=kM 



The zero lag of the LRACF gives the longitude-resolved vari- 
ance: 



erf = LRACF/o 

AT— 1 



7=0 



(5) 
(6) 



From this we computed the modulation index as m, = cr,-///,-. 
Since the zero lag of the LRACF is simply the sum of all el- 
ements of the corresponding fluctuation power spectrum, we 
may view the variance in a given longitude bin as having con- 
tributions from different frequencies (Parseval's theorem). We 
made use of this fact to compute modulation indices with the 
exclusion of the lowest frequency bin, since it often contained 
a strong component which, by virtue of its low frequency and 
longitude-independent modulation index, we attribute to inter- 
stellar scintillation. 

The 2DFS was derived through a similar division of the 
data sequence into M-pulse blocks. Since it integrates in longi- 
tude, we only include a designated on-pulse longitude interval 
in the analysis. The squared modulus of the two-dimensional 
DFT of such blocks gives the 2DF power spectrum, which was 



integrated over all blocks: 



N/M-l 



2DFS'™ = — V 

! nb u 

(k+l)M-l B-l 

V ^ s >_p-2*yR{iblB+jllM) 



l=kM b=0 



(7) 



where B is the number of longitude bins. When possible, to 
correct for the autocorrelation characteristics of the noise and 
interference, the raw 2DFS of an equal-sized "off-pulse" longi- 
tude interval (or the mean from several intervals) was computed 
and subtracted from the raw on-pulse spectrum. 

The 2DFS has two dimensions that are both dimensionless 
frequencies. In this work we refer to quantities in the pulse 
number- and longitude-associated axes as temporal and lon- 
gitudinal frequencies respectively, and plot the spectra with 
the longitudinal frequency axis as the abscissa. We usually 
use units of (subpulse) cycles per (spin) period (cpp) for 
these quantities, which in the case of longitudinal frequen- 
cies also corresponds to degrees of subpulse phase per de- 
gree of longitude. The 2DFS is periodic (due to aliasing with 
the finite sample interval) and symmetric about the diagonal 
(due to the real-valued input), obeying the relation 2DFS,y = 
2DFS( ± i +C B)(Tj+dM) where c and d are arbitrary integers, and can 
therefore be fully specified by computing a contiguous region 
of either 5/2+ 1 xM coefficients, oiBxM/2 + 1 coefficients. In 
this work we plot the region given by i e [-B/2: B/2)(i.e. lon- 
gitudinal frequencies between -18O°/A0 and 180°/A</>, where 
Acp is the longitude sampling interval) and j e [0 : M/2] (i.e. 
temporal frequencies between and 0.5 cpp), which generally 
avoids having responses "wrap" across borders in the planar 
projection, and directly gives information concerning the ap- 
parent sense and rate of any subpulse drift pattern. 

The two-dimensional DFT of the 2DFS gives a form of 
Two-Dimensional Autocorrelation Function (2DACF), a func- 
tion which, like the LRFS, has seen frequent use i n the past for 
drifting subpulse analysis (e.g. lTavlor et aljfl975h . Again, it is 
important to pad the data with zeroes (in both axes) to avoid 
computing a cyclic convolution, and to discard odd-numbered 
rows and columns of the resultant 2DFS to avoid misinter- 
pereting the significance of spectral features. The correction for 
edge normalization effects was performed through multiplica- 
tion by 1 -j/M; since the pulsar signal is not expected to extend 
beyond the longitude boundaries, additional normalization by 
1 - i/B is not appropriate. The result is: 



2DACF = 



1 



NB(l-j/M) 

N/M-l {k+\)M-\-j B-l-i 
Z Z Z S 'bI S (b+iKl+j)- 



(8) 



fe=0 l=kM b=0 



The zero pulse -number lag of th e 2DACF gives the single-pulse 
ACF used bv lJenet et all i boOll) : 



ACF, = 2DACF/0 

N-l B-l 



~ M ZZ^'^Uoj- 

7=0 b=Q 



(9) 
(10) 
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As an additional diagnostic of pulse-to-pulse correlations, 
we computed the LRCCF, which measures the covariance be- 
tween every pair of longitude bins as a function of pulse num- 
ber offset. We define the LRCCF as 

j N-k-l 

LRCCF,,* = — J] S' hj S' hU+k) . (11) 
j=o 

Our definition differs somewhat from that of Popov (1986), 
who normalized by l/cr^cr^ to obtain a true correlation rather 
than a covariance. We find the covariance preferable in the eval- 
uation of the significance of low S/N results since the noise 
levels are identical in all bins of each two-dimensional (11,1*2) 
map. The computation of the LRCCF could, in principle, be 
incorporated within the above scheme by integrating the con- 
jugate product of complex LRFS pairs and taking the DFT 
of the result, however the computational complexity is much 
greater than for the autocorrelation techniques above. Since we 
were only interested in the first ~ 10 (<§: M) lags we found 
it more efficient to compute directly from Eq. \H] We note 
however that the frequency domain counterpart (the Longitude- 
Resolved Fluctuation Cross-Power Spectrum) could be useful 
in simultaneous multi-frequency observations for examining 
the presence, strength and relative phase of periodic modula- 
tions between different longitudes and frequencies. 

Finally, although individual pulse detections were not the 
focus of this work, the assembled data set provides an opportu- 
nity for a survey of giant pulse behaviour in recycled pulsars. 
For this reason, we searched each observation for significant 
single pulses by convolving each pulse with boxcars of widths 
of 1, 2, 3, 4, 6, 8, 12 ...Wmax bins (where W max is approxi- 
mately the width of the on-pulse region of the average profile) 
and selecting those pulses containing convolved samples with 
a single-trial significance exceeding 6<x. We then computed the 
energy of the selected pulses as the sum of intensity values over 
an interval centered on the same longitude as the detection box- 
car but with twice the width, and expressed these as a fraction 
of the energy of the average profile. 

3. Results 

3.1. Non-Detections 

In the observations of PSR J0218+4232, J0034-0534, 
J0613-0200, J1643-1224 and B1937+21, with the exception 
of very low frequency fluctuation attributable to scintillation, 
we found no detectable modulation in any of the measured 
statistics. As with all pulsars observed, this included searches 
of the LRFS for transform lengths from 2 9 to 2 17 pulses (unless 
the number of pulses observed was less than this), from which 
one may derive limits on the presence of quasi-periodic mod- 
ulation. Considering a longitude-resolved variance computed 
over a given interval Av, the noise terms in the sum result in 
an approximately Gaussian distribution (central limit theorem) 
with a standard deviation proportional to Av 1 / 2 . For a given 
noise level estimated from the off-pulse parts of a LRFS, we 
take the square root of the 3-cr point of the noise in a given 
estimated variance (Eq.|6|l as the minimum detectable standard 
deviation. These values are given for each observation in Table 



^ The minimum detectable standard deviation for a given fre- 
quency integration is given by cr^Av/cpp) 1 ' 4 , which we have 
tabulated as cr m for the full modulation index (Av = 0.5 cpp) 
and for unresolved (i.e. nearly coherent) features in the longest 
transform, as <x c (Av = 1/M max cpp, M max = 2 17 unless noted). 
With the exception of PSR B 1937 +21, the four pulsars for 
which we have not detected modulation are those with the 
worst detection limits, and for this reason the results cannot 
be taken as evidence for low modulation indices. 

The search for bright individual pulses was successful in 
several pulsars (noted below), however with the exception of 
PSR B1937+21, no pulsar produced clearly detectable giant 
pulses in our observations. For PSR J0613-0200, one pulse 
barely exceeded the threshold, occurring within the (broad) on- 
pulse region, with an energy exceeding the mean by a factor 
of 20 and a width of one bin. Further studies should be made 
before attributing this detection to a giant pulse phenomemon. 
For all pulsars, the 6-cr energy limit for individual samples is 
presented in Table [2 as E m i n ; for wider pulses the limit scales 
as « 1/2 where n is the number of samples to add. 

3.2. PSRJ1012+5307 

The fluctuation spectra for PSR J1012+5307 (Fig.Q indicate 
the presence of quasi-periodic modulation in all three major 
profile components. Moreover, it is apparent in the 1380 MHz 
observation that the peak frequency of the modulation varies 
with longitude, showing different values in components la, lb 
and II. To better examine the differences between the compo- 
nents, we have plotted fluctuation spectra from components la, 
lb and lib in Fig.|2] From this figure it is apparent that although 
the overall depth of modulation in component la is higher than 
than that in lb, its quasi-periodic feature is broader in frequency 
space. We also note the presence of narrow features in the spec- 
tra of components lb and lib at frequencies close to twice the 
frequency of peak power. The significance of these features de- 
teriorates with a narrower smoothing function, indicating that 
their intrinsic width is close to that seen in Fig. [2] 

We detected 70 individual bright pulses in the 1380 MHz 
observation, the strongest of which had an energy five times 
the mean pulse energy for this obervation. The pulse profile 
formed by adding these pulses shows relatively more emission 
in components la and lib compared to lb and Ha, consistent 
with the higher modulation indices found in these regions. 

The two-dimensional fluctuation spectra were contami- 
nated by the presence of broad stripes aligned with the tem- 
poral frequency axis, which we were unable to subtract due 
to the lack of sufficiently wide off-pulse longitude intervals. 
Nevertheless, it is clear that the modulation power in the com- 
ponent I is divided between two main regions, one with zero 
frequency in the longitude axis, and one with a freqency of 
~ -10 cpp (Hig. [3j. The latter component indicates the pres- 
ence of quasi-periodic drifting subpulses, with a longitude 
separation of ~ 40° between subpulses. Since the compo- 
nent I is double-peaked, one might also consider the alterna- 
tive of longitude-stationary subpulse modulation with a sim- 
ple phase offset between the two halves of the profile com- 
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Pulse longitude (deg) Pulse longitude (deg) 

Fig. 1. Longitude-resolved modulation statistics for PSR J1012+5307 at 1380 MHz (left) and 840 MHz (right). Top panel: mean 
pulse profile (top, solid lines), standard deviation (unfilled circles) and modulation index (filled circles, right-side axis scales) for 
PSR J1012+5307 at 1380 MHz (left) and 840 MHz (right). The names used for components are indicated in the top left panel. 
The standard deviation and modulation index are subject to a 3-cr threshold in the power values from which they are derived, 
to ensure their significance. A strong component in the DC frequency bin of the LRFS is attributed to scintillation and was not 
included in the calculation of these statistics. Bottom panel: LRFS, shown as a greyscale computed with M — 512, along with 
contours of Gaussian-smoothed version of the same spectrum of FWHM 0° x 0.02 cpp (1380 MHz) and 20° x 0.02 cpp (840 
MHz). 



ponent. However, this can be rejected because the peak fre- 
quency obtained could not be more than 180° divided by the 
component separation of ~ 40°, i.e. ~ 5 cpp, where there 
is certainly no peak in the spectrum. The spectral component 
centered around zero longitudinal frequency indicates the pres- 
ence of longitude-stationary quasi-periodic modulation, while 
its relatively narrow extent suggests it extends across most of 
the component I. 

It is conceivable that the two spectral components in the 
2DFS are in fact associated with with same modulation feature, 
if the apparent temporal frequency offset is actually spurious. 
In this case the shifted, double-peaked nature of the spectrum 
in the longitudinal axis is indicative of linear drifting with a 



frequency of ~ 5 cpp and a phase step between opposite halves 
of the profile component fedwards et all2003l) . 

3.3. PSRJ1022+1001 

This pulsar showed non-periodic modulation at 328, 840 and 
1380 MHz (Fig.B). At 1380 MHz, the trailing component is 
seen to be very strongly modulated, a fact corroborated by the 
dominance of this component in the average profile derived 
from the ~ 10 4 individual pulses we detected. In this observa- 
tion, the brightest pulses had an energy ~ 6.5 times the mean, 
while the distribution of energy measured in a longitude in- 
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Fig. 2. Fluctuation spectra for PSR J1012+5307 at 1380 MHz, 
from peak bins in components la (solid), lb (dashed) and lib 
(dotted; Ha appeared similar and is omitted for clarity). The 
raw spectra were convolved with a Gaussian of FWHM 0.02 
cpp before plotting. 
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Fig. 3. 2DFS for component I of PSR J1012+5307 at 1380 
MHz, using M — 512 (greyscale) and smoothing by convo- 
lution with a Gaussian of FWHM x 0.02 cpp (contours). 

terval spanning the trailing component was consistent with an 
exponential or power-law. 

The fluctuation spectra of the bins near the two peaks of 
the average profile at 1380 MHz are shown in Fig. |5] While 
the leading component shows a flat fluctuation spectrum, the 
trailing component proceeds smoothly to a maximum at the 
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Fig. 4. Average profile, standard deviation and modulation in- 
dex for PSR J1022+1001 at 1380 (top), 840 (middle) and 328 
MHz (bottom). See caption for Fig.^ but here the modulation 
index is plotted with the same scale as the other quantities. 

Nyquist frequency. This corresponds in the autocorrelation 
function to a strong negative spike at a lag of one pulse. There 
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Fig. 5. Fluctuation spectra from peak longitude bins of the lead- 
ing (open circles, dashed line) and trailing (filled circles, solid 
line) components of PSR J1022+1001. The points show the 
measured spectrum (M = 512) in the peak bin, while the solid 
and broken lines show the spectra from the peak and surround- 
ing bins respectively, after smoothing by convolution with a 
Gaussian of FWHM 0.02 cpp. The low-frequency excess is due 
to very low freqency modulation (attributed to scintillation), 
broadened by the Gaussian smoothing kernel. 



is no evidence for any periodicity associated with this correla- 
tion. 
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3.4. PSR J1 51 8+4904 



Pulse longitude (deg) 

Fig. 6. Longitude-resolved statistics for PSR J1518+4904. See 
caption for Fig.^ 



We detected quasi-periodic modulation throughout the double- 
peaked main profile component of PSR J1518+4904 (Fig.|6j, 
with a frequency of approximately 0.38 cpp. In the peak bin 
this modulation contributes about half of the total modulation 
power, however this fraction decreases significantly toward the 
trailing part of the profile. The 2DFS (Fig. [7} shows a cor- 
responding component with its peak centered at a longitudi- 
nal frequency ~ -10 cpp, corresponding to a phase slope of 
about half a cycle across the pulse window. The 2DACF (Fig. 
[H} confirms this description, showing the characteristic repeat- 
ing pattern of drifting subpulses. The fact that for any given 
pulse number lag the 2DACF is singly-peaked confirms what 
one might expect with the relatively small phase slope, namely 
that usually only one subpulse is present as part of the quasi- 
periodic pattern. Fig. |9] shows this more clearly for the zero 
lag, and confirms that the subpulses are significantly narrower 
than the average profile. The drift pattern is also manifest in 
the LRCCF (Fig. HOi. from which it is clear that drifting of the 
modulation applies over the entire longitude window defined 
by the average profile. 



We detected 1650 individual pulses from this pulsar, the 
brightest of which had an energy 4.8 times the mean. 

3.5. PSR J1 71 3+0747 

We have detected modulation in observations at PSR 
J1713+0747 from 840 to 2240 MHz (Fig.QT)- The evolution 
in the pulse profile and modulation index is very minor over 
this range, however to aid in the assesment of the longitude- 
dependence of the modulation index under the low longitude 
resolution available, we have provided plots for all four fre- 
quencies. It is clear that the leading edge shows a higher mod- 
ulation index than other parts of the profile. A selection of 
one-dimensional fluctuation spectra are shown in Fig.^] Two 
broad maxima at frequencies of ~ 0.17 and ~ 0.35 cpp are 
clearly seen. 

We detected ~ 50 individual pulses at 1 190 MHz, with en- 
ergies up to 6 times the mean. 
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Fig. 7. 2DFS 
for Fig. 



Fluctuation frequency (cycles/period) 

for PSR J1518+4904 at 1355 MHz. See caption 




Pulse longitude lag (deg) 

Fig. 8. Lags 0-24 of 2DACF for PSR J 15 18+4904 at 1355 
MHz. Black and white values correspond to the most positive 
and negative values respectively, outside the zero lag. 



3.6. PSR J1 91 8-0642 

We detected modulation in the strongest two longitude bins in 
an observation of PSR J1918-0642 at 1380 MHz (Fig. [jj). 
There is some evidence for a concentration of power between 
0.25 - 0.5 cpp in the fluctuation spectrum (Fig. I14L however 
the limited sensitivity makes this result uncertain. 




10 15 20 

Pulse longitude lag (degrees) 



25 



30 



Fig. 9. Single-pulse ACF (solid line) and ACF of the average 
profile (dashed line) for PSR J1518+4904 at 1355 MHz. 

3.7. PSRJ2145-0750 

In observations at 860 and 1380 MHz, this pulsar showed 
strong intensity fluctuations, at some pulse longitudes exceed- 
ing 100% (Fig. [15}. Two 1380-MHz observations provided 
sufficient sensitivity to detect weak quasi-periodicities around 
~ 0.22 and 0.45 cpp, we plot the result from the better observa- 
tion in Fig. ^] To check for stationary or drifting modulation 
in the extended bridge and trailing component, we computed 
the the 2DFS of this longitude interval. No evidence for quasi- 
periodicity was found, but very broad features such as those 
seen in the peak bins cannot be ruled out. In the 1380 MHz 
observation we detected 490 individual pulses, the strongest of 
which had an energy 3.7 times the mean. The profile formed 
by adding these pulses was dominated by a narrow component 
around pulse longitude 215°, indicating that intensity enhance- 
ments in this component are not accompanied by enhancements 
at other longitudes. The less sensitive 840 MHz observation 
yielded similar results for the 1 1 detected pulses. 

4. Discussion 

The results of this study show that the emission from recycled 
pulsars strongly varies on timescales comparable to the rota- 
tion period. Of the seven pulsars for which we were sensitive 
to (aperiodic) modulation with m > 0.5, we detected such mod- 
ulation in all but one (PSR B 1937+21). Quasi-periodicities also 
seem to be a common feature of the fluctuation spectra, how- 
ever features are typically broad and often do not comprise 
the dominant contribution to the observed modulation index. 
Such a broad, weak feature has also been observed in PSR 
J0437-4715 llVivekanandetalll998h . 

For PSR J1012+5307 and PSR J1518+4904 the presence 
of offset components in the two-dimensional fluctuation spec- 
tra indicates that the phase of the quasi-periodic modulations 
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Fig. 11. Average profile, standard deviation and modulation index for PSR J1713+0747 at (a) 840 MHz, (b) 1 190 MHz, (c) 1700 
MHz and (d) 2240 MHz. See also caption for Fig.Q] 
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Fig. 12. Fluctuation spectra from significant bins of the LRFS, at 1 190 MHz (left) and 1700 MHz (right), for PSR J 17 13+0747. 
See also caption for Fig. [5] 
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Fig. 13. Average profile, standard deviation and modulation in- 
dex for PSR J1918-0642 at 1380 MHz. See caption for Fig. 
□ 



advances from the leading edge to the trailing edge of the main 
profile component. That is, the components of the emission that 
are periodically modulated appear at successively later longi- 
tudes. In addition to the results presented here, we checked 
other observations and found similar behaviour. These are the 
first ever detections of drifting subpulses in recycled pulsars. 
For the other pulsars the sensitivity to modulation at longitudes 




0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.4 0.45 0.5 
Fluctuation frequency (cycles/period) 

Fig. 14. Fluctuation spectrum from the peak longitude bin of 
PSR J1918-0642, at 1380 MHz. A Gaussian of FWHM 0.05 
cpp was used to produced the smoothed spectrum. (See caption 
for Fig. 0) 



away from the peak is probably too poor to have allowed de- 
tection of similar phenomen a, if present. For P S R J04 37-4715 
the sensitivity achieved by IVivekanand et al.l |l998) appears 
more than sufficient to have ascertained whether the modula- 
tion drifts, but unfortunately the two-dimensional autocorrela- 
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Fig. 15. Longitude-dependent statistics for PSR J2145-0750 at 860 MHz (left) and 1380 MHz (right). See caption for Fig.ffl 
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Fig. 16. Fluctuation spectrum from the peak longitude bin 
(solid line) and the bins immediately before (dashed) and af- 
ter (dotted) it, for PSR J2 145-0750, at 1380 MHz. See caption 
for Fig. 



tion function was computed without prior subtraction of the av- 
erage profile, causing it to be dominated by this component and 
potentially making any weak drift pattern invisible in the plot. 
Determining whether or not the subpulse modulation drifts 
woul d provide a good test of the model of ICril & Krawczvkl 
which predicts stationary modulation. 



Despite the very different physical conditions in the po- 
lar cap regions and magnetospheres of recycled pulsars ver- 
sus ordinary pulsars, the observed single-pulse phenomena are 
broadly consistent. Specificially, like ordin ary pulsars they tend 



Il97ll iBartel etal.ll 1980 


Weisbers et al. 1986S) and non- white 


fluctuation spectra (cf. [Taylor & Huguenin 197lUBacker et alJ 



distinguishes between modulation in so- 
called conal and core components of the emission, noting that 
conal components frequently show quasi-periodic modulation 
with periods of 2-10 pulses, which can sometimes be highly 
coherent, while core components tend to possess weak, broad 
low frequency excesses or peaks in their fluctuation spectra 5 . 
We clearly do not detect any features of high coherence in 
the recycled pulsars studied here, however it is difficult to say 
what proportion of normal pulsars show them, for studies have 
tended to favour pulsars known to possess them over those with 
broad components that are more difficult to characterize. None 
of the pulsars studied here appears to show a low-frequency 
excess (except for the very low frequency scintillation compo- 
nents), perhaps indicating that recycled pulsar emission is not 
related to core emission in slow pulsars. 

An apparent exception to the above generalizations is PSR 
B1937+21. Our non-detection (m < 0.26) is consistent with 
the pr eviously measured modulation index of 0.18 djenet et alJ 
2001 ) 6 . The reported constancy of the modulation index with 
longitude (except for the trailing edge, where giant pulses oc- 

5 iRankinl 1 1986t> draws attention to a number of reports of low- 
frequency quasi-periodicities in core-associated emission (P 3 ~ 20 
pulses), however examination of the source literature (or at least, those 
that published their spectra) indicates that with some exceptions, most 
are just as consistent with a low frequency excess as with quasi- 

periodicity_. 

2001) 



6 Since Ijenet et alJ y001) dealt with baseband samples, the ob 
served variance included the radiometer noise of the pulsar sig 
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cur) suggests scintillation as a possible origin for the modula- 
tion, and lJenet et al show that the level of modulation is 
consistent wit h the known prope rties of diffractive scintillation 
in this source. ljenet et al.ld200lh use this result and the similar- 
ity of the single-pulse ACF to the ACF of the average profile 
to argue that the intrinsic emission of PSR B1937+21 is very 
stab le, in contrast to all previously studied pulsars. However, 
IWeisberg et alJ ( fl9*8 6) find that pulsars they classified as "core" 
had low modulation indices, and indeed their result for PSR 
B2053+36 resembles that of PSR B1937+21 not only in the 
low value of modulation index but also in its constancy across 
most of the profile and in the shape of the profile itself. The 
issue could be resolved by examining the frequency structure 
of the modulation seen in these pulsars. 

In ordinary pulsars, the presence of drifting subpulses is 
usually associated with a viewing geometry that traces the 
emission along a q uasi-transverse path near the edge of the 
beam (e. g.lRank inll986l ). This is consiste nt with the picture put 
forth by Ruderma n & Sutherland! d 19751) that the pulsar beam 
consists of a circulating ring of subbeams corresponding to 
such a pattern of spark activity on the polar cap. The obser- 
vation that drifting subpulse patterns are coherent between dif- 
ferent pulse longitudes and between observations made simul- 
taneously at different radio frequencies (i. e. that the temporal 
period P3 is independent of these; see e.g. iTavlor et aljfl975l) 
lends strong support to this claim, as does the recent appar- 
ent detecti on of mo dulation at the cir culation period in PSR 
B0943+10 fcteshpande & RanlrinbOOll) . The presence of fluc- 
tuation components at different temporal frequencies ought, in 
such a model, to obey a mirror symmetry between pulse lon- 
gitude pairs that sample the same part of the rotating subbeam 
system. A clear example of this is PSR B1237+25, where the 
relative importance of the (at least) three modulation features 
shares the same approximate symmetry in longitude as the av- 
erage profile (Bac keJl973l) . 

How PSR J1012+5307 fits into this picture, if at all, is un- 
clear. Taken at face value, the variation of peak fluctuation fre- 
quency with pulse longitude implies that components la, lb 
and II all correspond to different (nested) rotating subbeam 
cones. The much higher modulation index of component II ver- 
sus component I would support this assertion. In this case one 
might expect that the emission from 290 to 340° corresponds to 
the complementary components required to preserve symmetry 
mentioned above and ought to show a progression from low to 
high fluctuation frequencies with longitude (i.e. in reverse to 
that seen from la to lb to II). 

The detections of individual bright pulses made here are 
also broadly consistent with the behaviour of ordinary pul- 
sars, the pulse energy distributions of which typically extend 
to a factor of a few higher than the mean (e.g. Ritchings 
Il976l). A simi lar result for PSR J0437-4715 was presented 
bv lJenet etall dl998t) . With the exception of PSR B1937+21, 



nal itself. The observed modulation index (for which they reported 
1.032 ± 0.001 and 1.034 ± 0.003 in the main pulse and interpulse 
respectively) is then given by (1 + 2m 2 ) 1 ' 2 . In the use of integrated 
samples this additional variance is reduced by the reciprocal of the 
time-bandwidth product and is typically neglected. 



for which giant pulses are already known to occur, and PSR 
J0613-0200, for which the result of a single detection needs 
further investigation, we did not detect any giant pulses. If 
any of the pulsars observed had produced unresolved giant 
pulses with the same energy probability distrib ution as that re - 
ported for PSR B 1937 +21 at 430 MHz JCognard et aDl 19961) . 
we should have detected in t he worst case (PSR J02 18+4232) 
~ 90 of them. The results of Kinkhabwal a & ThorsetJ §000) 
appear to indicate that for PSR B 1937 +21, giant pulses are 
less frequently observed at higher frequencies, although they 
do not provide a fit to the p robability distribution. Our rough 
estimate based on Fig. 8 of Kinkhab wala & Thorsettl (|2000), 
given a mean flux density of 10.4 mJy (from the spectral fit of 
lFosteretalJfl99lh and a ratio of 0.68:0.32 in energy between 
the main pulse and interpulse (measured from our observa- 
tions) is that giant pulses of energies exceeding a given multiple 
of the mean energy occur of the order of 30 times less often at 
1420 MHz than at 430 MHz. Giant pulses with a similar proba- 
bility distribution would have been detected if present in any of 
the other pulsars observed. This result is particularly interesting 
for the first three pulsars listed in Tabled since among field re- 
cycled pulsars the inferred magnetic field strengths at the light 
cylinder for these pulsars are exceeded only by those of PSRs 
B1937+21 and B1957+20, making the m prim e candidates for 
giant pulses according to lCognard et alJ Jl996l) . 

Other than the detections of giant pulses in PSR B1937+21 
and PSR B18 12-24, we note that there exist in the litera- 
ture largel y unexamined reports of strong pulses from PSR 
B1534+12 JSallmen & Backerll995tlSalrmerJl998l) . The pulse 
energy distribution clearly extends to much higher energies 
than observed here for other recycle d pulsars (exc ept PSR 
B1937+21), with the version given by Sallmen ( 1998) depict- 
ing the detection of pulses up to 17 times the mean energy, from 
a sample of just 12000 pulses. One wonders whether, were a 
data set of a number of pulses comparable to those used in stud- 
ies of PSR B 1937 +21, PSR B 182 1-24 and the Crab pulsar an- 
alyzed (~ 10 7 ), sufficiently bright pulses would be received for 
them to be considered "giant". Such a result would be very in- 
teresting, since the inferred magnetic field strength at the light 
cylinder for PSR B1534+12 is low compared to most recycled 
and young pulsars (~ 1.6 x 10 3 G, versus ~ 10 6 G for PSR 
B1937 +21 and ~ 9 x 10 s G for the Crab pulsar). The sugges- 
tion of Cognard et al. ll 19961) that this parameter is an indica- 
tor of giant pulse activity could be strongly tested by studying 
strong pulses from PSR B1534+12. 

Given the lack of a complete theory of pulsar emission, 
there are few predictions that can be tested by the new obser- 
vations presented here. A recent attempt at a systematic pre- 
scription for the properties of emission as a function of pulsar 
paramete rs for both normal and millisecond pulsars is that of 
iGil & Sen dvk (2000). They postulate that the polar gap region 
is populated by a system of sparks that drift about the center 
as a set of concentric "rings". Assuming that the spark width 
and separation are both equal to the gap height, they define a 
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"complexity parameter" (a), equal to the ratio of spark width to 
the width of the polar cap, which they calculate as 



a = 5 



10" 



2/7 



(n) 



-9/14 



(12) 



where P is the pulsar spin period. Since the space between 
sparks is equal to their width, the total number of sparks across 
the cap is ~ a, and the number of nested cones of emission 
(excluding the centre spark) is n (a - l)/2. They show that 
the complexity parameter is correlated with profile morphol- 
ogy, with core single profiles having the highest complexity 
parameters, followed by triple profiles and then multiple and 
conal single profiles. This leads them to suggest that core single 
profiles, and the central components of triple profiles, are the 
result of emission from a large number of narrow cones, broad- 
ened by the finite beam- width of the elementary emission (2/y, 
where y is the Lorentz factor of the emitting particles). This ex- 
plains why their profile morphology is simple, and makes the 
prediction that the modulation index of such pulsars should be 
small. 

In Figure El we plot lines of constant complexity param- 
eter along with the period and intrinsic period derivative for 
46 field pulsars presumed to be recycled (P < 0.1 s, P < 
10~ 15 ) . We used the publ i shed p roper motions o flToscano et alJ 
i 19991) and lLange et"al] ll200lh . plus a new measurement for 
PSR J15 18+4904 Qx = 7.5 ± 0.7 mas yr _1 ; Stappers un- 
published ), in combination w ith the Galactic potential model 
of Kuiiken & Gil mord Jl989l) to correct the observed period 
derivatives f or the effects of acceler ation with respect to the 
line of sight toamour & Tavl or 1991). These points are shown 
as solid circles, while uncorrected points are not filled. The 
sources for which the modulation properties are now known 
are labelled. For most sources the model predicts 1-3 cones 
of emission, implying complex profiles with significant mod- 
ulation, as observed. The exception to this statement is PSR 
B 1937 +21, with a large complexity parameter of a a 23. This 
pulsar is also exceptional for its low modulation index and sim- 
ple profile morphology, just as expected for a large number of 
unresolved cones. We have already drawn an analogy between 
PSR B1937+21 and PSR B2053+36, and it is interesting to 
note that the complexit y parameter of the lat ter is also quite 
high: a 10. Indeed, Weisber g et alJ II1986I) found that not 
only PSR B2053+36, but all of the core-dominated pulsars in 
their sample have low modulatio n indices, in excellent agree- 
ment with the assertion, made by Gil & Sendvk (2000) on the 
basis of the tendency for core pulsars to have large complex- 
ity parameters, that core components are the sum of numerous 
1 /y-smeared components. 

We have shown that the detection and characterization of 
pulse-to-pulse intensity fluctuations in millisecond pulsars is 
possible using existing instruments. The acquistion of similar 
information about southern recycled pulsars and those accessi- 
ble from Arecibo should expand the available sample and make 
stronger generalizations about the population possible. With 
the possible future availability of very large radio telescopes 
such as LOFAR and SKA both the number of recycled pulsars 
within the reach of study, and the depth of detail to which they 
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Fig. 17. Period vs Period derivative for 46 recycled pulsars in 
the Galactic disk. The solid points have been corrected for ac- 
celeration effects, the unfilled points have not. The dotted lines 
correspo nd to lines of e qual value in the complexity parameter 
of lGil & Sendvkl foOOd) . a = 20, 10, 5, 2.5 (top to bottom). 



could be studied will be greatly enhanced and hopefully pro- 
vide considerable constraining information on recycled pulsar 
polar caps and magnetospheres. 
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